Simulating local measurements on a quantum many body system with stochastic 

matrix product states 
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We demonstrate how to simulate both discrete and continuous stochastic evolution of a quantum 
many body system subject to measurements using matrix product states. A particular, but generally 
applicable, measurement model is analyzed and a simple representation in terms of matrix product 
operators is found. The technique is exemplified by numerical simulations of the anti-ferromagnetic 
Heisenberg spin-chain model subject to various instances of the measurement model. In particular 
we focus on local measurements with small support and non-local measurements which induces long 
range correlations. 
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I. INTRODUCTION 

With the experimental realization of ultra-cold atoms 
and optical lattice-systems experimentalists begin to 
probe various many-body models accurately, leading to 
theoretical and accurate experimental studies of highly 
non-trivial physical ideas, such as topological states of 
matter, frustrated systems and phase transition phenom- 
ena. In these systems experimentalists are not only mea- 
suring macroscopic observables, but are also beginning 
to measure local observables [l|, H, Q and correlation- 
functions JiJ. 

Understanding quantum many-body systems is a chal- 
lenge due to the large dimensionality of the many-body 
Hilbert space and the associated complexity of the states. 
Not only can it be difficult to describe a generic many- 
body state, but finding the physically relevant states (e.g. 
diagonalizing the system Hamiltonian) , and quantifying 
their physical properties can be extremely difficult. 

Recent progress in the understanding of entanglement 
and the complexity of quantum many body systems have 
led to several simulation techniques for strongly corre- 
lated systems with local interactions based on matrix 
product states (MPS) for a one-dimensional lattice sys- 
tem, projected entangled pair states (PEPS) and vari- 
ants hereof in higher dimensions [1, @ 0, B H, OH EH, 
which enable us to calculate ground states, perform time- 
evolution of states and calculate expectation-values of 
many interesting physical operators accurately and effi- 
ciently. 

In this paper we will use the matrix product state 
techniques to simulate local and non-local measurements 
on quantum many-body systems. Measurements lead 
to interesting conditioned dynamics and provide alter- 
native routes to entanglement generation [l2l . Il3j and to 
quantum computing [14j . Quantifying interaction and 
measurement-induced dynamics occurring on the same 
time-scale is, however, highly non-trivial. 

In section|TI]we review the formalism and highlight the 
basic features of matrix product states. In section ITTTl we 
will briefly review the quantum theory of measurements 



Figure 1: (Color online) Schematic of the measurement setup. 
A two-level ancilla interacts for a time r with a specific site 
or a collection of sites in the lattice, after which a projective 
measurement is performed on the ancilla. 



and introduce the measurement model studied in this 
paper. In section IIVI we describe how to simulate both 
discrete and continuous measurements on matrix product 
states, and in section [V] we show some example simula- 
tions. Finally we conclude and discuss further work in 
section I VII 



II. 



MATRIX PRODUCT STATES AND 
OPERATORS 



We wish to study the effect of measurements on many- 
body systems but, as mentioned in the introduction, solv- 
ing many-body problems can be extremely difficult, and 
with the added complication of non-equilibrium stochas- 
tic behavior the problem is not getting any easier. We 
will describe how to simulate a measurement scheme in 
terms of matrix product states (MPS) [1, EJ or projected 
entangled pair states 10] (PEPS) (the methods presented 
here can be extended very easily from the latter to the 
former). It will turn out that local measurements and 
certain global measurements have a natural description 
in terms of matrix product states, and since measure- 
ments are intimately connected to entanglement, this de- 
scription will also highlight aspects of the entanglement 
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properties of matrix product states. 

First we will briefly review the basics of matrix product 
states. Any use of this technique begins with a factor- 
ization of the full many-body Hilbert space into elemen- 
tary constituents of dimension d. This factorization is 
done such that the states of interest only contain a lim- 
ited amount of entanglement between its subsystems (in 
the sense of small Schmidt-number between any bipartite 
cut). For spins on a periodic one-dimensional lattice of 
length L it is natural to factor the Hilbert-space into a 
product of single sites. The MPS ansatz is a parametriza- 
tion of the expansion coefficients in the product basis 

L 

|V> = Y. Tr (tW'^tPI'" 2 . . . T^^) (g) \a t ) , (1) 

ai...cr L i=l 

where T^ ,ai are d DxD matrices (or a dx D x D tensor), 
where D is called the virtual dimension and |(Tj) denote 
the d single-site basis states. The trace operation in ^ 
is carried out over the virtual dimensions and results in 
coefficients of the different product states — the larger the 
virtual dimension D, the more entanglement is supported 
by the ansatz. 

If one wishes to consider a system with open boundary 
conditions one can take the first and final matrices 
and T^l to have dimensions 1 x D and Dxl respectively. 
The essential feature of this parametrization is that any 
Schmidt-decomposition of the above state will contain at 
most D terms, thus limiting the entanglement entropy 
between any bipartite splitting in a systematic way. Re- 
ciprocally, it can also be shown that any state where all 
bipartite splittings contain at most D terms are of this 
form Q. 

Single site reduced density-matrices for states of this 
form can be calculated in a clever way [l5| using 
only 0(Ld 2 D 3 ) operations. Similarly, few-site density- 
matrices can be calculated efficiently. Thus, expectation 
values of local and correlation-observables can be calcu- 
lated efficiently for matrix product states. 

Using variational methods it is possible not only to 
find good approximations for ground-states of nearest- 
neighbor Hamiltonians, but also Greens-functions, ther- 
mal states and time-evolution of these states, etc. [8Hl5j. 
All these variational methods essentially rely on the ma- 
trix element (x|A|-0) between any two matrix-product 
states and an operator A to be calculated efficiently, 
which is often possible for physically relevant operators, 
like tensor-products of local operators. 

A special set of efficiently contractible operators are 
the matrix product operators (MPO), which are con- 
structed in the same way as matrix product states. An 
MPO is parametrized as 

L 

A= Tr(M [1] ' jl ...Af [i]ji ) 0I 3i , (2) 

jl,—,jL i=1 

where M''l Ji are DxD matrices and the N = d 2 op- 
erators Xj constitute a basis for the operators on a sin- 
gle site. If we choose the matrix-basis as the Xj (i.e. 



| cr) (ct'I), the matrices M' l ' Ji can be thought of as matri- 
ces of operators. 

If A is of this form, applying A to in 
([T]) amounts to the update rule ^(na')(0p') = 

T^pM^p (aAX^), where {oca') is to be un- 
derstood as a combined index as in the Kronecker matrix- 
product Cnty/jn — AijBki- The virtual dimension of 
A \ ip) will then be the product of the virtual dimensions 
of A and In practice however, it is not desirable and 
often not necessary to increase the virtual dimension as 
one can use a variational method 0] to find a matrix 
product state \x) of a given virtual dimension that mini- 
mizes || \x) — A \tf>)\\ ■ In, e.g., time-evolution, where A is 
an approximation to the unitary time-evolution operator, 
one is interested in keeping a fixed virtual dimension, and 
the error associated with this truncation is usually negli- 
gible for short times. As shown in [lj| it is also possible 
to perform the time evolution in the Heisenberg picture 
using matrix product operators. 

The above techniques can be generalized in terms of 
projected entangled pair states (PEPS) [lOj], where the 
state is written in terms of a general tensor-network 
[ItJ instead of just nearest neighbor contractions on a 
one-dimensional lattice. For these kinds of states it 
is possible to formulate variational methods in much 
the same way as for MPS, but the numerical stability 
and time-complexity of the methods can depend drasti- 
cally on the topology of the graph. For an open one- 
dimensional lattice with nearest-neighbor interactions 
calculating ground states and performing time-evolution 
scales as 0(LD 3 ), whereas, for a one-dimensional lat- 
tice with periodic boundary-conditions, the calculations 
scales as 0(LD 5 ) due to the loop in the factorization- 
graph [l8j . 

The topology of the factorization-graph is usually de- 
termined by the interactions in the system, and may be 
a chain, a loop, a grid, a tree or some variation of these 
depending on the system. Since non-local measurements 
can induce entanglement directly between the measured 
subsystems due to the obtained information, the topol- 
ogy can also depend on the type of measurements being 
performed on the system. 



III. THE QUANTUM THEORY OF 
MEASUREMENTS 

If we wish to simulate the stochastic evolution of a 
quantum system subject to measurements we can use sev- 
eral, closely related, techniques depending on the type of 
measurement being performed. 

When a quantum system is being measured, the quan- 
tum state \tp) changes because the measurement alters 
the observers knowledge of the system. Formally this 
back-action is described by the action of a set of opera- 
tors where the state conditioned on the measure- 
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ment outcome fi is fi M \ip) /y (^fi^O^I?/;). In general, 
such a set of operators is a valid measurement scheme 
if fij,f2p = ^ can be shown [l9| that the ac- 
tion of any such set of operators can be implemented by 
coupling the measured system unitarily to an auxiliary 
system and performing a projective measurement on the 
auxiliary system. In particular, consider a two-level an- 
cilla interacting with the probed system. If a z acts on 
the ancilla space and A is an observable of the probed 
system with spectral resolution A = J2 a aP a {P a is the 
projector onto the eigenspace associated with the eigen- 
value a), then an interaction on the form 

V = ga z A, (3) 

applied for a time r will result in a time-evolution de- 
scribed by 

U = exp(-W) = v ( 2( ^ a i z ) ® p a (4) 

a 

where T> is the spin 1/2 rotation-operator on the ancilla- 
spin system and = gr. In words, depending on the 
value of the system observable A, the effective spin an- 
cilla will be rotated an angle 20a around the z-axis. If 
we initialize the ancilla spin in state \+ x ) pointing in the 
x-direction and subsequently measure the y-component, 
we thus gain partial information about the operator A. 
Performing this measurement with the outcome fi = ±1 
for the cr ^-readout, corresponds to the action of the mea- 
surement operators 

^ = {(f, y \ ® 1)17(1+,.) ® 1) = i (e-^ A + ifie^ A ) (5) 

up to an arbitrary phase. Note that since we are coupling 
the system to a qubit-ancilla we can maximally gain one 
classical bit of information from the system for each mea- 
surement. Controlling cf> by varying the interaction time 
r or strength g we can vary the strength of the measure- 
ment. For (f> — no measurement is performed. Con- 
tinuous measurements with infinitesimal changes of \ip) 
correspond to the limit of frequent applications of weak 
measurements, <j) — > 0, while for ~ 1 the measurement 
back-action becomes large. 

We can use this formalism to propagate a system be- 
ing measured with some rate k by simply evolving the 
state \ip) unitarily with system Hamiltonian H, for a time 
St = 1/k, and then apply the measurement operator Q^: 
\ip) — > il M \ip) I \/p((a), where we assume the measure- 
ment is much faster than St and fi is chosen randomly 
according to the distribution p(/i) = (^\Qlfl fJ ,\'ip). This 
is then repeated until the desired final time. 

If the system is continuously monitored, we model the 
measurement process as a limit of the above procedure 
in the following sense: When all the fJ^ are infinitesi- 
mally close to the identity, and the measurements are 
performed with a large rate compared to the system dy- 
namics, the accumulated effect of measurements in any 



given time-interval (small on the timescale of unitary 
evolution) amounts to iV„ applications of f2^ where 
are binomially distributed random numbers. Performing 
the limiting procedure results in a stochastic differential 
equation. 

To be precise, if <f> <C 1 and 1/k ^ t sys then the accu- 
mulated effect of N = kAT measurements in a time AT 
is given by 

Q(N+, N-) oc 1 + 4>A{N+ - N-) 

+ ^<f> 2 A 2 ((N + - N_ f -2N). ((>1 

In the limit N ^> 1 we can apply the central-limit theo- 
rem to obtain N+ - N- « 2kAT0 (A) + y/nAW, where 
AW is a normally distributed stochastic variable with 
zero mean and variance AT. The accumulated effect of 
the measurements can then be written as 

n(N+, N-) cx 1 + 2AT0 2 k (A) A - -AT^nA 2 

2 ! ( 7 ) 

+AW^/ r ^A, 

Including the Hamiltonian evolution of the system and 
state-normalization, then in the limit AT — > dt we get 
the stochastic differential equation 

d |</>) = -iHdt \ip) + dt- (2 (A) A- A 2 - (A) 2 ^j \ip) 

+Vf{A-{A))dW\i>), 

(8) 

where T = (f> 2 K is the measurement strength. 

IV. STOCHASTIC PROPAGATION OF MATRIX 
PRODUCT STATES 

If we consider measurements on a system described 
by an MPS then we proceed as described in section [TTT1 
Propagate \ip) unitarily for a time St and then apply f2 M 
with probability p(fi) — (ip IW^D.^ip), If is efficiently 
contractible we can use the same variational principle 
as used for time-evolution, to find the matrix-product 
state of a given virtual dimension which best approximate 
^). 

In particular, if can be written as a matrix- 
product operator the matrix-element for any 
MPS ansatz \x) can be calculated. The probability 
p(fi) — flpl'ijj) can also be calculated efficiently 
since fij^f^ is also a matrix-product operator, although 
possibly of a higher virtual dimension. In particular if 
the system is measured at a single site, fi^ and f^fi^ 
are just tensor products of a number of identities and a 
single-site operator. In this case the application of f2 M to 
the MPS is of course trivial just as the application of any 
local operator V is trivial for an MPS: Simply update 

T Hf a ' = E< T a/' fal^K). where V acts on site i. 
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We thus find that many-body systems with local in- 
teractions combined with single-site measurements are 
easy to simulate using matrix product states in one di- 
mension or PEPS in higher dimensions. If, however, the 
measurement extends across multiple sites, e.g. if our 
probe is not absolutely confined to a single site, then we 
need to decompose f2 M into an MPO and use the vari- 
ational methods in order to avoid growth of the virtual 
dimension of the state. 

We will now consider the measurement discussed in 
section Iffll described by ©, © and ©, but with A a 
sum of local operators, i.e. A = J2i9jOj- This can arise 
either as a combined simultaneous interaction, but also as 
result of a sequential interaction of an ancilla with multi- 
ple sites. As in equation (O the measurement operators 
for this measurement is given by 




1(1 ( 



J4>9jOj 



(9) 



(10) 



We see that applying this f2 M to a general MPS re- 
sults in a superposition of two MPS each a copy of the 
original state but multiplied with the product-operators 
®j exp(±i<figjOj). One should think that such a super- 
position state would not be well represented by an MPS, 
but if we increase the virtual dimension of the state from 
D to 2D we can in fact represent the state exactly as well 
as any superposition of two MPSs [2(| • 

Indeed, fi M can be written as an MPO with tensors 



MW = 
AfW = 

= 



giOi) ii^exp(i03iOi)) 

A) \ 

exp(i<j>giOi) J 

( exp(-icj)g L L )\ 
\ exp(i(fig L L ) J ' 



(§exp(-i<; 
/ exp(-i<& 
I 



(11) 



1. e., of virtual dimension 2. Note that it is the use of a 
2-dimensiona! ancilla, that results in virtual dimension 

2. Had we used an n-dimensional ancilla, generically we 
would need a virtual dimension of n for the MPO. For 
the calculation of the branching probabilities p(/i), we 
note that also Qtfi M can be represented by an MPO of 
virtual dimension 3, since fttft M = i(l+2//i(exp(2# J 4)- 
exp(-2i</»A))). 

As mentioned above, the variational method can now 
be applied to find a matrix-product state with a spec- 
ified virtual dimension which is closest to fl^ \ip). The 
truncation-error for the posterior state will naturally de- 
pend not only on how close to the identity f2 M is, but also 
on the state being measured. As an example, an initially 
uncorrelated spin-chain where all the spins point in the 
x-direction, will get long-range correlations if we, e.g., 
measure the z-component of the total spin to 0. If how- 
ever, a measurement of the rr-component is performed, 
the spin-chain will remain in a product state. 



In the above measurement scheme the rate k is the in- 
verse of the simulation time-step and the measurement 
strength <f> can be adjusted. If we wish to model contin- 
uous measurements, the measurement rate and strength 
can no longer be chosen independently and the accumu- 
lated effect of many measurements for each simulation 
time-step has to be taken into account. One approach 
would be to implement the limit discussed in section IIIII 
by simply choosing smaller time-steps and scaling <f> ac- 
cordingly. But then a larger number of truncations will 
be performed, so even if one is able to apply the mea- 
surement exactly the regular time-evolution may lead to 
an increased error. 

If we seek to approximate the continuous measurement 
regime, the time-evolution is given by ([5]), which is not 
only non- linear in \ip) via the terms (A) A and (A) , but 
also contains terms proportional to A 2 , which cannot 
readily be cast into a matrix-product form. In each time- 
step, however, we can calculate (A) and pick a random 
AW to select the relevant il-operator. To apply this op- 
erator to \ip) we need to evaluate (x|f2|-0) as described in 
section|nl but even though O is not on a clear MPO-form 
one can just evaluate the terms separately, provided A 2 
is not too pathological. This amounts to the stochastic 
Euler method which has a global error of 0{At 1 / 2 ) (2~H ]. 

Alternatively, and in general, if the measurement ex- 
tends over more than one site, it is of course always pos- 
sible to combine those subsystems into a single Hilbert 
space and then apply the measurement operator directly 
on that space. If the measurement has sufficiently small 
support, this will not ruin the power of MPS, since it 
only treats a small part of the system exactly. 



V. NUMERICAL EXAMPLES 

We have applied the above techniques to a spin- 1/2 
chain with an anti-ferromagnetic Heisenberg Hamiltonian 



H 



(12) 



where the sum is over nearest neighbors. We first find 
an MPS for the ground state, and with this as our initial 
state we begin to probe the system with the measurement 
scheme outlined above. 



A. Local measurements 

The interesting regime in this numerical study is when 
the measurement- and Hamiltonian dynamics are of com- 
parable importance, but let us first imagine a projective 
measurement of the z-component spin on one of the end- 
points of the chain. Since the reduced density-matrix for 
each site in the anti-ferromagnetic ground-state are com- 
pletely mixed (the Hamiltonian is rotational invariant), 
the expectation value of the local spin is zero and hence 
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Figure 2: (Color online) Examples of posterior states after 
different types of measurement on the spin-chain ground state. 
In both plots <j> = 0.3-7r/4 using the measurement operators 
(0. The plotted values are {a\) (red) (cr*) (blue) and (p\o l z ) 
for the ground state (solid black) for states subject to f2 + (A = 
al) in[2]i and Q+(A = al)Q,+ (A = cr%) in[2p. 

the projective measurement will have probability 1/2 for 
both \+ z ) and |— z ) . The spin is, however, correlated with 
its neighbors, and the spin-expectation-value of the con- 
ditioned state (o"z) cond is directly related to the ground 
state correlation function. 

Rewriting the correlation-function between two arbi- 
trary spins in terms of the projector onto eigenspaces of 
the z'th spin-z, P+/_i wc obtain 

K4) = \ {(no (p-O) = = - (p-<) , 

since the ground state is invariant under a z — > — a 3 z for 
all j. But (P±a J z ) is proportional to the conditioned ex- 
pectation value of a{. 

To get a feel for what a single measurement (close to 
the projective case) can do to a spin-chain ground state 
see figure [2] Figure [2^ shows data for the state 0+ 
where fi M is given by ((5]) with cf> — 0.37r/4 and A = a\. 
We see, as noted above, the spins align (on average) with 
the z-z correlation function and nothing happens to the 
Cy-average. In addition, {p\a\) is less modulated than 
in the initial ground state due to the partial projection. 
In figure [2)3 we show the effects of a joint measurement, 
where we couple two ancillas to the system — one to a\ 
and one to er^. Notice how the z and y-averages follow 
their respective correlation-functions from each end of 
the spin-chain. 

If the system- and measurement-dynamics occur on 
the same timescale the dynamics becomes far more com- 
plex. We have simulated repeated measurements over 
time for a spin-chain with 60 spins with varying measure- 
ment strengths (f>, as shown in figure [3] where time-series 
for (erg) is shown. For weak measurements (black solid 
and dashed) the measurements are not strong enough 
to project the spin to a er 2 -eigenstate, since the spin is 
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tJ 

Figure 3: (Color online) Trajectories for single-site measure- 
ments with varying strengths showing both the quantum 
zeno- and weak-limit of the measurements. In all simulations 
k — 100J and the values of <j> are shown in the legend. 



strongly driven by its neighbors. For strong measure- 
ments (red dashed) the interaction is not strong enough 
to drive the measured spin away from its measured value 
quickly enough, and we observe a quantum zeno ef- 
fect on the first spin; effectively pinning the spin to 
a random, but definite, direction. In the intermediate 
regime (red solid) the measured spin exhibits oscillations 
with comparable signatures from both measurement- and 
interaction-dynamics. 

In figure[2]the single-site scheme has been simulated (as 
in figured) and the expectation- values of <r z for the first 
three spins are shown as function of time. As noted above 
the measurements induce oscillations in the average-spin 
values, and in the example shown, the measurements are 
strong enough to project the first spin onto an eigenstate 
for a long interval of time. Note also that the second and 
third spins reflect anti-correlation and correlation with 
the first spin respectively, both in time-intervals, when 
the spin is almost in a <r z -eigenstate, and when {a\) is 
close to zero. 

Notice, that the anti-correlations persist even after the 
measured spin has been almost completely projected. 
This can be understood, if we consider the weak mea- 
surements as a stochastic perturbation: In order to excite 
the high-energy non-anti correlated eigenstates of H the 
stochastic perturbation needs to have significant support 
at frequencies of comparable magnitude. Since this par- 
ticular measurement is weak it does not provide enough 
energy to excite these states. 

Figure [5] shows simulations, where a\ measurements 
are performed on site 1 and <jy at the end L = 60 as 
in figure ^Bp. Here we also see the characteristic oscilla- 
tions of measurements competing with interactions. In 
this case the measurements appear to be completely inde- 
pendent. In the very long time-limit, one might expect to 
see temporal correlations arise due to propagation-effects 
along the chain. 
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Figure 4: Time-series of (al) (solid), (<r%) (dashed) and (erf) (dash-dotted) for measurements using {(Sj with A = a\ for . 
and k = 100J. 
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Figure 5: (Color online) Time-series for measurements using © with both A = a\ and A = for <j> = 0.05 and k = 100 J. 
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Figure 6: (Color online) To deal with non-local measure- 
ments, e.g., by sequential system-ancilla interactions (left), 
in the PEPS-formalism it is necessary to add an addi- 
tional entanglement-bond in order to provide an efficient 
parametrization of the many-body state (right). 



B. Non-local measurements 

In figure [5] we consider a process where we measure two 
sites far from each other at the same time. This might 
be accomplished by an ancilla first interacting with site 
i, then with site j and then being projectively measured. 

This can be simulated in two ways within the theory 
of MPS and PEPS. One way is to keep the usual matrix 
product states and simply apply © in the form of MPO 
as in (jlip where — for k ^ The result of such 
a simulation can be seen in figure where we have sim- 
ulated a measurement of the sum a], A + a^ 1 in a lattice 
of 30 spins. In the simulation we see many of the same 
features as in figures [4] and [5l where the spin-expectation 
values are projected onto eigenstates of the observable 
A. In this case the observable A also has eigenstates 
with eigenvalue 0, and in the simulations (figure [7]) we 



see a projection onto this eigenvalue for t > 4. By closer 
inspection of the simulation data it was found that in the 
periods where the measured spins have average value 0, 
the two-site density matrix is a slightly asymmetric sta- 
tistical mixture of the \S = 0, m = 0) and \S = 1, m = Q) 
state with a purity of about 50%. 

For sites not too far apart this works well, but 
truncation-error will quickly become a problem due to 
the emerging entanglement between site i and j. If 
the sites are far apart, however, one might introduce a 
PEPS-graph as shown to the right in figure [6l where an 
extra entanglement bond has been introduced between 
the two distant sites. This topology constitutes a new 
parametrization of our states, but it is also a good illus- 
tration of the physical effects of the measurement: We 
are now, effectively dealing with a kind of closed bound- 
ary conditions mediated by measurements. This entails 
an increase in computational complexity as mentioned 
above, but it is still possible to use many of the same nu- 
merical tricks and techniques as used in open- and closed 
boundary conditions in usual MPS-simulations. 



VI. CONCLUSION AND OUTLOOK 

In summary we have shown how to simulate measure- 
ments on a quantum many-body system described by a 
MPS or PEPS formalism. In particular, a natural class of 
measurement-schemes can be represented simply by ma- 
trix product operators and their PEPS-equivalents. We 
have illustrated the use of these techniques on the anti- 
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Figure 7: Time-series for a measurement using ([5]) for A = cr 14 + a^ 1 , <j> — 0.1 and k — 100 J. 



ferromagnetic Heisenberg spin-chain ground state and in- 
vestigated the dynamics resulting from particular weak 
measurement schemes. 

There are a number of natural applications and exten- 
sions of this work. Cold atoms in optical lattices con- 
stitute a very attractive model of many-body dynamics, 
where both direct optical imaging by a high aperture lens 
and by the transmission properties of an optical cavity 
enclosing part of the atomic ensemble [H, [HI, US HH HE] 
is possible. High resolution achieved through non-linear 



atomic response [26J, |2jj as well as localized ionization 
signals, due to the impact of a scanning electron beam 
[l|, may be modelled by our approach. 

This will allow studies of the interplay between mea- 
surement induced and interaction induced localization 
phenomena in such models. In a future perspective, 
closed feedback-loops on quantum many body systems 
may constitute a promising application to perform more 
general quantum state engineering and possibly to con- 
trol phase transitions in many-body systems. 



[1] P. Wurtz, T. Langen, T. Gericke, A. Koglbauer, and 
H. Ott, Phys. Rev. Lett. 103, 080404 (2009). 

[2] W. S. Bakr, J. I. Gillen, A. Peng, S. Foelling, and 
M. Greiner, arXiv:0908.0174 [cond-mat] (2009). 

[3] I. B. Mekhov and H. Ritsch, Phys. Rev. Lett. 102, 020403 
(2009). 

[4] T. Rom, T. Best, D. van Oosten, U. Schneider, S. Foiling, 

B. Paredes, and I. Bloch, Nature 444, 733 (2006). 
[5] Y. Y. Shi, L. M. Duan, and G. Vidal, Phys. Rev. A 74, 

022320 (2006). 
[6] G. Vidal, Phys. Rev. Lett. 101, 110501 (2008). 
[7] F. Verstraete, D. Porras, and J. I. Cirac, Phys. Rev. Lett. 

93, 227205 (2004). 
[8] G. Vidal, Phys. Rev. Lett. 93, 040502 (2004). 
[9] F. Verstraete, J. J. Garcia-Ripoll, and J. I. Cirac, Phys. 

Rev. Lett. 93, 207204 (2004). 
[10] F. Verstraete and J. I. Cirac, arXiv:0407066 [cond-mat] 

(2004). 

[11] A. J. Daley, C. Kollath, U. Schollwock, and G. Vidal, 
Journal of Statistical Mechanics: Theory and Experi- 
ment 2004, P04005 (2004). 

[12] D. N. Matsukevich, P. Maunz, D. L. Moehring, S. Olm- 
schenk, and C. Monroe, Phys. Rev. Lett. 100, 150404 
(2008). 

[13] A. S. S0rensen and K. M0lmer, Phys. Rev. Lett. 90, 
127903 (2003). 

[14] R. Raussendorf and H. J. Briegel, Phys. Rev. Lett. 86, 
5188 (2001). 



[15] F. Verstraete, V. Murg, and J. I. Cirac, Adv. Phys. 57, 
143 (2008). 

[16] M. J. Hartmann, J. Prior, S. R. Clark, and M. B. Plenio, 

Phys. Rev. Lett. 102, 057202 (2009). 
[17] J. I. Cirac and F. Verstraete, arXiv:0910.1130 [cond-mat] 

(2009). 

[18] I. L. Markov and Y. Shi, SIAM Journal on Computing 
38, 963 (2008). 

[19] M. A. Nielsen and I. L. Chuang, Quantum computation 
and quantum information (Cambridge University Press, 
2000). 

[20] I. P. McCulloch, Journal of Statistical Mechanics: The- 
ory and Experiment 2007, P10014 (2007). 

[21] P. E. Kloeden and E. Platen, Numerical solution of 
stochastic differential equations (Springer, 1995). 

[22] M. Karski, L. Forster, J. M. Choi, W. Alt, A. Widera, 
and D. Meschede, Phys. Rev. Lett. 102, 053001 (2009). 

[23] A. Takamizawa, T. Steinmetz, R. Delhuille, T. W. 
Hansen, and J. Reichel, Optics Express 14, 10976 (2006). 

[24] I. B. Mekhov and H. Ritsch, Phys. Rev. A 80, 013604 
(2009). 

[25] I. B. Mekhov and H. Ritsch (2009), arXiv:0911.0389 
[quant-ph] . 

[26] D. D. Yavuz and N. A. Proite, Phys. Rev. A 76, 

041802(R) (2007). 
[27] A. V. Gorshkov, L. Jiang, M. Greiner, P. Zoller, and 

M. D. Lukin, Phys. Rev. Lett. 100, 093005 (2008). 



